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Abstract. We study morphologies of thin-hlm diblock copolymers between two flat and parallel walls. 
The study is restricted to the weak segregation regime below the order-disorder transition temperature. 
The deviation from perfect lamellar shape is calculated for phases which are perpendicular and parallel 
to the walls. We examine the undulations of the inter material dividing surface and its angle with the 
walls, and find that the deviation from its unperturbed position can be much larger than in the strong 
segregation case. Evaluating the weak segregation stability of the lamellar phases, it is shown that a 
surface interaction, which is quadratic in the monomer concentration, favors the perpendicular lamellar 
phase. In particular, the degeneracy between perpendicular and unfrustrated parallel lamellar phases for 
walls without a preferential adsorption is removed. 

PACS. 61.25.Hq Macromolecular and polymer solutions; polymer melts; swelling. - 61.41. +e Polymers, 
elastomers, and plastics. - 68.55.-a Thin films structure and morphology. 

1 Introduction field theory, the bulk phase diagram is governed by two 

parameters: / = Na/N, the fraction of the A-block in a 

, n , ,T->/~n->\ i ri chain of polymerization index N — Na + -/Vr, and yN, 

Diblock copolymers (BCPJ are made up ol two chemically 1 J ' A ' 

i _li i i i . i, mi nnTi where y is the Flory parameter measuring the interaction 

distinct chains covalcntly bonded together. IhcBCPsys- ° 

c nr , , j , , , between the two species, and is inversely proportional to 

tern forms sell assembled structures with length scales m r ' j r * 

the nanometer to micrometer range. On the level of mean- ^ e * em P era t urc 
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For temperatures above the order-disorder transition Previous analytic theories p^,p3|,p4[, while providingvalu- 
(ODT) temperature the system is in the disordered phase, able qualitative results, have been sensitive to the specific 
As the temperature is lowered, symmetric BCP melts (/ = choice of phenomenological coefficients, and this sensitiv- 
vj) undergo a weak first order transition to a lamellar phase ity leads to marked inaccuracies of the order parameter as 
at ^ > Xc- As the degree of block asymmetry / is in- compared to Monte-Carlo simulations [f22| . 
creased, |/ — || > 0, other phases of hexagonal and cubic In the present work we complement the numerical stud- 
spatial symmetries become stable ies by introducing an alternative analytical method. In 

m, . , r . , n , . r n/-(Ti 1 1 n particular, we derive the deviation of the perpendicular 

Ihe mterfacial behavior of BCP melts has been the 

subiect of exoerimental § U ITol Fill and theoretical O HI and P arallel lamellae from their bulk shape. In Sec. | we 

HHHHHIIIIig investigations. In the former case introducc a modcl free ener gy and derive the underlying 



, , , , n • j 1 , 1 r,r<r, j u equations. In Sec. M the shape of confined lamellae is in- 

the substrate is typically spm-coated by the BCP, and sub- L-l 

, n 1 n 1 ii t , lv vestigated and found to be, in general, very different from 

sequently analyzed by small angle neutron and X-ray scat- a 

n It ,, nl the bulk shape. The energy of the perpendicular and par- 

termg or neutron reflectivity measurements. It the walls OJ ^ 

, . ., , e ,. , n , . r allel lamellae as a function of surface separation 2L as well 

are neutral, i.e., without preferential adsorption to one 01 

, , , , . cl c -1 nl J., , n , as the stability diagram is discussed in Sec. 0. 

the two blocks, thin films ol lamellar diblock copolymers jo u 

maintain their bulk periodicity do by aligning perpendic- 
ular to the confining walls. Such long range ordering can 
be transferred by various techniques to a surface, creating 
a template useful in nanolithography [ pT| . In cases where 
the walls prefer one of the two blocks, the lamellae can 



2 Model 



Close to the phase transition point (ODT) between the 
disordered and lamellar phases, the free energy of symmet- 
ric BCP melt is well described by the following Ginzburg- 

reduce the interfacial interactions by aligning parallel to , . j-i r— 1 1— i r— 1 r— 1 r— 1 

Landau expansion ||7|, |l4|lq ,|25|,[26|,p7||: 

the walls, and change the lamellar periodicity from its bulk 
value do. Which of the two phases prevails (parallel or per- 
pendicular) depends on the distance 2L between the two 
walls (film thickness), the strength of the wall interactions The bulk free energy T h (in units of k B T) is given as a 



^ 2 + i/,(V^ + , o 2 0) 2 + ^ 4 }d 3 r (1) 



as well as the degree of segregation N\- functional of the local order parameter </>(r) = cj>A (r) - /, 

which is the deviation of the A monomer concentration 



Numerical calculations of confined BCP have been per- 
formed using self-consistent field theory [jl5|[l(|] and Monte- 
Carlo simulations jl6],[2^]. Using these techniques, order 
parameter profiles and phase diagrams have been obtained. / = 1/2 



from its average value. The parameters above are given 
by 
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<7o — 1.9b /R g ; r = 2pN (x c — x) (2) const., used throughout this paper, no additional surface 



With a monomer size a, the gyration radius for Gaussian 
chains is Rg ~ ^Na 2 , and p = 1/Na 3 . Other parameters 
in Eq. (|l|) are 



length scales are introduced (see Refs. |lj,|l8|,[19[|29| where 
er(r s ) varies on the walls). The strength of wall interaction 
is given by two parameters: a + = a(y = L) and a~ = 
a(y = —L). Symmetric (er + = <t~) and asymmetric (<t + = 
Xc — 10.49 /N ; h = 3pc 2 R 2 /2q 2 (3) — er _ ) walls will be considered as special cases. 

The deviation of the order parameter, 0(r), from its 

The dimensionless parameters u/p and c are of order 

bulk value 0b (r) is denoted by 5<j> 

unity. The reduced temperature r ~ (xc ~ x) is positive in 

the disordered phase, where <j)(r) = 0. Close to the ODT S(f>(r) — <fr(r) — </>&(r) (5) 



This deviation contains the effect of the walls. The free 



the bulk system is described by two length scales: the first 

is the periodicity of lamellar modulations do — 27r/go, an d 

, . . , . . . , , ,, u . energy T — Th + T s is then expanded to second order 

the second is the correlation length ~ (r/h) 1 , charac- 



around its bulk value, T = T\4>h\ + A!F[S(f) 7 



+ -(r + -M 2 ) (St) 2 + -h {V 2 5cp + q$64>) j d 3 r 
+ f [a5<P + t s (2<j) b 5dp + 5q> 2 )] d 2 r s (6) 



terizing the decay of surface induced modulations. This 

length diverges at the ODT, y = y c . . _ f \ r/ , a. , 1 ,o , 9 „ 9 , , r . 

8 S AF = / | [(r + hqt)(j> b + -ucpl + hq 2 V 2 <b b }50 

The interaction free energy of a BCP melt with the , , ^ 

confining wall (in units of ksT) can be written as a sum 
of two terms 

■p s = J [ o -(r s )0(r s ) + T s 2 (r s )] d 2 r s (4) In the next section we investigate the parallel and per- 

pendicular lamellar phases denoted as L|| and respec- 
tively, and choose the appropriate forms for their bulk 
phase (fit,. The free energy, Eq. (Q), is then minimized with 
respect to the correction field 8<j) and yields the BCP pro- 
file. 



where {r s } denotes the wall position. The first term is 
linear in the order parameter, and expresses preferential 
adsorption: a positive &(r s ) induces a negative 4>(v s ) (pref- 
erence to the B monomers). The second (quadratic) term 
allows surface deviation of the Flory parameter x from its 
bulk value. A positive t s means that the surface has an 

ordering temperature lower than the bulk one @. 3 Order parameter profiles 

In the following we consider a thin film in which the 
melt is confined by two flat and parallel walls at y = ±L. 
Interactions between the wall and the melt are assumed 
to be short-range, and for homogeneous walls, cr(r s ) = 



The cases of parallel Lu and perpendicular L± phases are 
now considered separately. 
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3.1 Perpendicular lamellar phase: L± where g± = g{±L) and w± = w(±L). 

The amplitude function g(y) results from the surface 



Up to this point 5<p and <pb were not specified. For films 
below the ODT (t < 0), the perpendicular bulk phase L± 
has the bulk periodicity do = 2ix /qq. Its order parameter 
is given in the single mode approximation (close to the 
ODT) by |0) 

<pb(r) = 4> q cos(q Q x) (7) 

the amplitude <j) q = (—St/u) 1 / 2 is obtained from a varia- 
tional principle of the bulk free energy. 

The order parameter for the perpendicular lamellae is 



modification of the Flory parameter, t s 0, and it van- 
ishes if t s vanishes. This can be seen by noting that if 
t s = then the minimum of the integral in Eq. (|l^) is 
obtained for g(y) = (recalling that r < 0). There is no 
coupling between w(y) and g(y), since the free energy is 
expanded to second order in 8<p, and the mixed terms are 
of higher order. The function w(y) minimizes AF W sub- 
ject to the condition that J w(y)dy is fixed. Using A as the 
Lagrange multiplier, it satisfies an ordinary fourth order 
4>j_(r) = 4>t (r) + 6<f>(r) differential equation 



6<f>(r) =w(y)+g(y)cos(q x) (8) ^4 _ l_) w {y) + 2qlw"{y) + w""(y) - ~ = (12) 

where for the correction field 5<j> we use the single mode Similarly, the equation for g(y) is 
ansatz. If additional modes are included in the bulk order 2r 
parameter, Eq. (0), such modes should also be included in 
Eq. (§). 

For the above choice of <pb [Eq. (0)], it is now possible 
to perform the x and z integration explicitly, retaining 
only the y dependency in Eq. (Jfy. The free energy per 
unit area can be written as 



T .9 + .9"" = (13) 

A Lagrange multiplier is not needed here because 
J g{y) coa(qox) d 3 r = 0. These equations are linear in w(y) 
and g(y) since the free energy, Eq. (^|), is expanded to 
second order around cf>b- The four boundary conditions for 
g(y) are 



AF ± = AF g + AF W (9) 



where 



and 



2T s <f> q + 2t s9 ± t hg"'(±L) = (14) 

g"{±L) = (15) 

and for w{y) 

cr ± + 2t s w± T q 2 hw'{±L) qp hw"'(±L) = (16) 

+ T s <j> g (g+ + g^) + -r s (g 2 _ + g 2 + ) (10) 

2 qlw±+w"{±L) = (17) 

J ^w(y)dy = (18) 

AF ™ = J | -\ Tw2 + \ h (Vo™ + w") 2 | dy Equation @ expresses the condition that the total A/B 

(w 2 _+w\) (11) fract ion is conserved, J S(f>(r)d 3 r = 0. 



Y. Tsori and D. Andelman: Diblock copolymer thin films: Parallel and perpendicular ... 5 
All coefficients in Eqs. (|l^) and ( |T3| ) are constants, and surface (IMDS) (denned as the surface where <j>(x, y) = 0) 



therefore the solutions g(y) and w(y) have the form: is perpendicular to the walls. This result is similar to the 

w(y) = A w e~ k ™y + B w e k -y + A* w e~ k ^ + B* w e k ^ 



one obtained in Ref. |16 (their Fig. 3), using different 



+ const. (19) 
](y) = A g e- k »y + B g e k "y + A* e - k > + B*e k > (20) 



methods. In Fig. 1(b) we show the response fields g{y) 

and w(y) in 56 = w(y) + g(y) cos(qox). It is advantageous 

for the lamellae to reduce their amplitude close to the 

u j a -lu i j- r ii • i , wall, hence, in our convention, a positive r. > induces 

where O denotes the complex conjugate of the variable ' ' ' v 8 

O. The complex amplitudes A w , B w , A g and B g and the a ne S ative & = °)" The am P litude of sinusoidal modula- 
constant term in Eq. © are determined from the bound- tions in v) = w ^ + + ^v)) is therefore 

ary conditions diminished from it unperturbed value <j) q . In the absence 

mi i i j j i of surface fields, o , the w part of 86 vanishes, w(y) = 0. 

I he complex wavevectors k w and k g are given by ' ^ r ' yy/ 

Figure 2 is similar to Fig. 1, but the symmetric walls 

*&=-«g + (r//0 1/a (21) 

(cr + = cr ) are chosen here to favor the B monomers (in 

k 2 g = (2r/h) 1 / 2 (22) 

dark), which partially wet them. As a result, the A/B 
In the vicinity of the ODT, r w 0, and the real and imag- IMDS bends and intersects with the walls at an angle 
inary parts of k w = k' w + ik'^ are given approximately which is different than 90°. The preferred adsorption is 
by also seen in Fig. 2(b), where w(y) is negative at the walls, 

w± < 0. 

k' w *^(N X -N Xc ) 1/2 (23) 

/ at at \ The copolymer contour lines are defined by the relation 

1 // ~ I 1 2 ^' X * * Xc \ fr,A\ 

K « go I 1 - a — a J (24) 

\ / ^>(r) = 0(,(r) + r ) = c, where c is a constant. Clearly, 

where a = 2q 2 /(1.95\/3c) follows from Eqs. (g) and (||). the inter-material dividing surface (IMDS) is just the spe- 

The period of modulations 2-Kjk'^ tends to 27r/go, and the cial case with c — 0. For bulk lamellar phase the IMDS 

decay length of these modulations £ w — l/k' w diverges as are just parallel planar surfaces (lines in two dimensions), 

£ w ~ (Nx — Nxc)^ 1 ^ 2 in the limit r — ► [[l7|[ll|. but for lamellae confined in thin films the shape of these 

A contour plot of the order parameter <j>(x,y) = lines is more complicated. Figure 3(a) shows the IMDS 

6i,(x,y) + 56(x,y) is shown in Fig. 1(a), for inter-plate lines for symmetric walls, both favoring the B monomer. 

separation 2L = 3c?o- The two walls at y = ±L are neutral, As expected, the contact area of the B domains with the 

o~^~ = 0, but the surface Flory parameter deviates from its wall is increased, and the IMDS lines are curved appropri- 

bulk value, r s > 0. Note that the interfacial width broad- ately. A different behavior is seen in Fig. 3(b) (asymmetric 

ens close to the wall, but the A/B inter-material dividing walls) where the curving of the IMDS lines is opposite at 
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the two surfaces because of the opposite wall interaction. \N = 11, and the contour lines show stronger undula- 

The deviation from a perfect lamellar shape is seen as the tions. Close to the ODT, the length scale associated with 

IMDS undulates. the undulation periodicity is 27r/fc^, ps do [see Eq. (p4j)], 

In general, contour lines do not run perpendicular to but it may get much smaller as the temperature is reduced, 

the wall but rather form an angle different than 90° with X Xc- The second length scale in the system, 2ir/k' w , 

the surface. On contour lines having <p — const., x and y characterizes the decay of modulations, and it diverges at 

are related by the ODT. 

c~w(y) 

cosq x = - — ; — TT ( 25 ) 

09 + g{y) 

3.2 Parallel lamellar phase: L 

Figure 4 is a schematic presentation of the IMDS. The 

dotted vertical line shows the unperturbed location of the The alternative case of lamellar order occurs when the 

A/B IMDS. At the y = -L wall, the deviation Ax of the lamellae are parallel to the walls, and the A/B profiles de- 

IMDS from this line (see Fig. 4) is pend only on the distance from the walls, 0(r) = <j>(y). In 

. the strong stretching approximation [ pT| , |20[ |, the lamellae 

— = — arccos(—— ) - - (26) 

a o Z7r <Pq t .9- 4 are allowed to stretch or compress in order to vary their 

The departure from the flat interface can be quite large, natural periodicity d according to the constraint inter- 

for example, in Fig. 3(a) it is Ax/d « 0.1. We define 9 plate separation 2L. The system can have n or n ± 1/2 

as the angle at which the IMDS line y(x) joins the wall lamellae between the walls, where n is the closest integer 

at y = -L. For neutral walls, Ax = and 9 = 90°. From to 2L/d . This strong stretching calculation motivates our 

Eq. (^|) it follows that zeroth order approximation to the L\\ phase, 

dy q Q sirL{q Ax){ct) q +g_) 2 <t>b{p) = ±4>q c °s[q(y + L)} (28) 



tan 9 = ^ = - *" " ; ' (27) 

dx w'(-L)(cj> q +g_)+w-g'(-L) 



Using the same parameters as in Fig. 3(a), we find that 
9 re 80°. 

As the ODT is approached from below (r < 0), the 
lamellae can be deformed more easily. The energetic cost 
of lamellae bending and compression is reduced, and the 
IMDS departs appreciably from its flat shape. The effect 
of temperature is clearly seen in Fig. 5, where in (a) the 
IMDS is plotted for Flory parameter x^ = 12, while in 
(b) the temperature is higher and closer to the ODT, 0||( r ) — 'Pbi^) + Sip(r) 



Using n from above, the wavenumber is q — nit / L or 
(n + 1/2)tt/L. The lamellae are stretched if q < qo and 
compressed if q > qo. The ± sign of the profile is deter- 
mined by the wall interactions. 

The bulk approximation for the profile, Eq. (p8[), serves 
as a starting point. However, the correction field, 8<f>, has 
an important contribution in the weak segregation. The 
order parameter for the parallel phase L|| is 



as: 
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S(p(r) — w(y) (29) recurrence relation converges only for specific values (eigen- 

values) of k. If k is a valid eigenvalue, then so are k*, 
— k and — k* . These four eigenvalues correspond to the 
four independent solutions of the fourth order differential 
equation, Eq. (|3ll). 



The free energy (per unit area) has y-dependent terms 
only. Expanded to second order in w(y), it can be written 



AF« 



^+2^+ \h (q 2 w + w'f 



d/y 



A useful approximation to the free energy Eq. ( |30| ) is 
+ a W- + 2T s (pb(—L)w- +t s w_ obtained by replacing </3| by its average, smeared value 

+ <j + w + + 2t s 4>] 3 (L)w + + r s tui (30) {4>1)- This is equivalent to replace the potential term 



Similar to the treatment of the perpendicular phase in 



Sec. 3.1, this free energy is minimized to yield a linear 



— (r/h) cos 2qy in Eq. (31) by its zero average. The gov- 
erning equation for the correction field w(y) is then given 
by a linear differential equation with constant coefficients 



differential equation, but with y-dependent coefficients: 



w""(y) + 2q 2 w"(y) w""(y) + 2q*w" (y) + (q* - r/h)w(y) - \/h = (36) 



9o - T ~ T ™s(2qy) w(y) - = (31) 



T T 



h h 



A 



h 



Under this approximation, the form of w{y) in the paral- 

The conditions imposed on w(y) are: , i T , ■ n • • ,i j- i T 

1 Kt " lei Lii phase is the same as it is m the perpendicular L± 

+ „ , / , rs « phase, Eq. dl3), only the boundary conditions are diffcr- 

a ± + 2r s (f) b {±L) + 2t s w± 1 1 H J J 



T qlhw'{±L) T hw"'{±L) = (32) 



cut. 

Order parameter profiles are presented in Fig. 6. The 
dotted line is <pb as obtained by the bulk approximation, 

qlw{±L) + w"{±L) = (33) 

l Eq. (p8|), and the solid line is the full profile, = <frb(y) + 

/ w(y)dy = (34) 

J—L &4>{y)- I n Fig. 6(a) the interfacial interactions are the same 

where as before A is the Lagrange multiplier and the last on both walls, <t+ = er~, inducing a symmetric lamellar 
equation expresses the conservation of the relative A/B ordering. The difference between the two curves is the cor- 
concentration in the film. The homogeneous solution of rection field S<f)(y), favoring adsorption of the B monomers 
Eq. @ has the Bloch (Floquet) form {<t> < 0) at the two walls. In 6(b) the film is asymmetric 

with <j + = —<j~, and adsorption of the A monomers at 

w{y) = e- kv ^a n e 2mqv + e- k " v ^a* n e- 2mq v (35) 

n n the y = —L wall is enhanced. 

A recurrence relation between the coefficients {a n } is ob- 
tained by substituting Eq. (|35| ) in Eq. ([jl"|). However, the 
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4 Free energy and Stability diagram proximation calculation. However, the difference between 



Once the order parameter profiles for the parallel and per- 
pendicular lamellar phases are calculated, the correspond- 
ing free energies can be evaluated by substituting the or- 
der parameter profiles in Eqs. ( |30| ) and (^), respectively. 
The reference free energy F° [(pi,] is calculated by the bulk 
approximation. For the parallel lamellae it is given by sub- 
stituting the profile, Eq. (p8|), directly into Eq. (|l]), 



the two curves tends to zero as 2L/cIq — ► oo, because the 
surface induced modulations have finite range. For both 
choices of cr, the L± free energy is significantly lowered 
from its bulk approximation value. Note that the bulk ap- 
proximation curves are similar to the curves obtained in 
the strong stretching approximation. 



Restricting ourselves to Ln and L± lamellar phases, 
the stability diagram is constructed as a function of two 
system parameters: the inter- wall separation 2L and the 
where the ±er terms result from the choice of the bulk surface preference a ± . In the first stage, we ignore the 



2L 

± a-(/> g ± <x + <f> q + 2r s <^\ (37) 



order parameter (j> b . For the perpendicular phase Lj_, sub- correction presented above, and use F9[<j> b ] and Ff[(j> b ] as 
stituting the profile , Eq. (0), results in: given by the bulk approximation calculation. The stability 



pirn 



1 ±2 , u ±4, 

-rd>„ H 6„ 

4 Vq 64 Vq 

2Lt 2 /-8t 



2L + t s $\ 



(38) 



diagram in Fig. 8 is calculated for walls having a fixed ra- 
tio of surface interaction er + = — 2a~ . Parallel lamellae at 
2L/do = n, for integer n, have symmetric ordering, while 



As a function of inter-plate separation 2L, the total antisymmetric ordering occurs for 2L/d^ = n + 1/2. The 

free energy F + 5<p] has oscillations, as depicted in difference in the diagram is caused by the choice of er's. In 

Fig. 7 for symmetric film, o~ = a + . The free energies of 8(a) the surface Flory parameter is the same as the bulk 

the perpendicular and parallel lamellar phases (solid line one, t s = 0. For neutral walls, a ± = 0, the perpendicular 

and thick solid line, Eqs. (||) and (|30|), are lower than the lamellae (in dark) are stable. A degeneracy between L± 

bulk ones F° [4>b] and F± [4>b] (dotted and dashed lines). In and L\\ phases occurs for 2L/do = n or n + \ (n = 0, 

Fig. 7(a) the wall interactions are — 0AhqQ<f) q , and the 1, 2 ...), where the parallel lamellae are not frustrated 

free energy of the Lh phase is slightly reduced from the and q = q . The parallel lamellae (in light color) are pre- 

bulk approximation value. Additional minimum develops ferred if the surface interaction is strong enough to over- 

at 2L « 1.5c£o- The L± free energy has a marked minimum come the lamellar stretching or compression. The use of 

for 2L < do see inset. In Fig. 7(b) the interfacial in- a weak segregation bulk approximation agrees with pre- 

teractions are smaller, = Q.2hq^(j) q , and in this case viously obtained stability diagrams in intermediate and 

the Lii free energy is notably lowered from the bulk ap- strong segregations |lq,^q|. 
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In Fig. 8(b) we present the bulk approximation, but 5 Conclusions 

now the surface Flory parameter is changed, r s = Q.lhq^ > 
0. The L|| phase is pushed upward and the diagram is dif- 
ferent. Symmetric phases [2L ndo] are pushed more 
than the asymmetric ones [2L « (n+ ^)do) because of our 
choice of surface fields er's. In the bulk approximation, the 
free energy of unfrustrated parallel lamellae [q — qo in Eq. 
([37])] is higher than that of the perpendicular lamellae [Eq. 
(|38|)] if the walls are neutral. As a result, the L± morphol- 
ogy is favored for all separations 2L, and the degeneracy 



is removed [15,16|. Clearly, a surface segregation temper- 
ature different than the bulk one (r s ^ 0) can account for 
the experimental lack of this degeneracy According 
to the same reasoning, if r s < then the parallel phases 
are preferred on the expense of the L± phases, and, in 
particular, for wall separations 2L sa ndo- Note that the 
last sentence agrees with the different dependence on r s 
seen in Eqs. ([37]) and (|38|). For parallel lamellae, the sur- 
face term is 2r s 4>j 1 , while for perpendicular lamellae it is 



only r s 



/,2 



We have used a Ginzburg-Landau free energy to study 
analytically the thin-film ordering of diblock copolymers 
(BCP) in the weak segregation regime. The two homo- 
geneous confining walls are assumed to have short-range 
interactions with the BCP blocks. The free energy is ex- 
panded to second order around the appropriately chosen 
bulk phase, and the correction field 8<f> is obtained. The 
use of such free energy formulation is advantageous be- 
cause it offers simple analytical results and complements 
numerical studies. However, our mean-field approach is 
limited to a region of temperatures in the vicinity of the 
ODT, but not too close to it, where critical fluctuations 
are known to be important ||. Very close to the ODT, the 
response field 8<j) diverges. However, if the surface Flory 
parameter is different from the bulk one, t s > 0, the sur- 
face has a lower ordering temperature than the bulk, and 
this divergence is absent p2] |. 

For confined parallel L|| and perpendicular L± phases, 
the correction field 8(f) adds an enrichment layer of the 



Figure 9 shows the stability diagram, where in (a) we preferred compon ent, with thickness £ ~ l/k' diverging 
use the bulk approximation for symmetric a~ = a+ film, at thc QDT Thig thickness is obt ained as a special case 
and in (b) we used thc full, and correct, order parameter for patterned walls ( in homogeneous a) studied by us be- 
<f> = b + 6<p. The parallel ordering is then always symmet- fore Effects of finite chain lengthj howeverj preclude the 
ric. A general feature of this diagram is that the L ± phase divergence of this thickness. In the L x phase, an increase 
is more stable relative to the L {] for larger surface fields of the surface fidds ff ± increases the correction field 5<j>, 
for 2L > d . The figure also demonstrates the qualitative and induces a paralM lameUar ordering untilj eventually, 
agreement with the bulk approximation. there ig nQ dear distinction between L {] and L ± . 

In general, the IMDS lines are bent and deviate from 
their flat shape in bulk lamellar system. Previous works 
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used a phenomenological model valid in the strong seg- In experiment with neutral walls, perpendicular lamcl- 

regation regime, and obtained a linear equation for the lae are always favored over unfrustrated parallel lamel- 

deviation of the IMDS. The resulting order parameter ex- lae (of period do) pi] ], in contrast to the common strong 

pressions for the confined phases are crude, when com- stretching prediction |0,n5lEn]. We first compute the 



pared to Monte-Carlo simulations [[^|23j,|24j. In the weak bulk stability diagram and find it similar to previous inter- 
segregation presented here, the order parameter itself is mediate and strong segregation calculations |15lnfl]. We 



linearized. Using the expressions given above for <pb (r) and then show that proper account of the surface change of 

8<f>(r), one can deduce the shape of an arbitrary equi-</> line the Flory parameter (r s > 0) can explain the experimen- 

given by </>(r) = c. We give expressions for the angle of the tal findings, and significantly change the stability diagram 

IMDS with the confining walls, and the deviation of the [compare Fig. 8(a) to 8(b)]. Thus, perpendicular lamellae 

IMDS from its flat shape. This deviation, characterized are expected to have the lowest free energy at all separa- 

by decaying oscillations, can be quite large and can even tions 2L, as in Fig. 8(b). We point out that if the surface 

reach 20% — 30% of the lamellar width <1q. We note that in ordering temperature is higher than the bulk ODT tem- 

an experimental setup whose target is to produce perpen- perature, i.e. t s < 0, the Lu phase may become stable even 

dicular lamellae, system parameters should be tuned in for neutral walls at 2L = ndo- However, this is yet to be 

order to keep the lamellae as flat and parallel as possible, confirmed experimentally. 



The free energy as a function of wall separation 2L 
is different from the bulk approximation. The free energy 
of the L± phase is lower than the one obtained the bulk 
approximation, as is seen in Fig. 7(a). The curve has de- 
caying oscillations and tends to a constant when 2L 3> e?o- 
The correction to the L» free energy has similar undula- 
tory character and under different conditions its effect can 
be large, as in Fig. 7(b). The pressure, —dF/dy, is differ- 
ent than what is expected from the bulk approximation, 
since additional maxima and minima are present in the 
free energy. Our bulk approximation yields order param- 
eter and energy profiles which are the same (apart from 
numerical values) as those obtained by the strong strech- Partial support from the U.S. -Israel Binational Founda- 



The stability diagram in this paper is similar to the 
diagram in [fl5|| . For symmetric walls, a + = <r~ , the L± is 
found to be stable for larger a fields than the bulk approx- 
imation predicts, while for 2L « do it is stable for smaller 
a fields. 

One possible way to refine the calculation presented 
here is to use a more accurate ansatz for the bulk order 
parameter <\>\, (r). Such ansatz will include more g-modes 
or an amplitude other than <j) q , further lowering the free 
energy. 
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7641 (1998). 

32. For a BCP at y > confined by a surface at y — with a 
uniform interaction parameter a, the correction field w(y) is 
w(y) = A w e~ kwV +c.c. The solution of Eq. © subject to 
boundary conditions analogous to Eqs. (|l6|), ( |l7j ) and ([L8|) 
then gives the amplitude 

A w — —ak w (4iT g k'^ + 2ig 2 (— r /ft)" 1 / 2 ) -1 , which diverges at 
the ODT unless r s ^ 0. 





Fig. 1. Contour plot of the perpendicular lamellar phases be- 
tween two homogeneous walls. The A monomers are shown in 
light shades while the B ones are dark. In (a) the two walls 
at y = ±L = ±1.5do are neutral, a ± = 0. Part (b) shows 
the correction fields w(y) (solid line) and g(y) (dashed line) 
in S(f>(r) = w(y) + g(y)cos(qox). A surface Flory parameter 
which is different from the bulk value, t 3 > 0, causes sur- 
face deviations of the lamellar structure from its bulk shape, 
even for neutral walls. The Flory parameter is X-^ = 10.8, 
t b — O.lhqQ. In this and subsequent figures we use u/p = c — 1 
and N = 1000. 
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(a) 




-1 1 



Y/d Q 



Fig. 2. Same as in Fig. 1, but here the two walls favor the B 
monomers, <r ± = 0.2hq?)(f> q > 0. Monomers are rearranged near 
the walls and the A/B inter-material dividing surface (1MDS) 
is curved (see also Fig. 3) . The preferential walls induce parallel 
ordering, as w(y) in (b). The length scale of modulation in 
(a) is determined by the functions w(y) and g(y) in (b) [Eqs. 
© and (§)]. 



(a) 




Fig. 3. Parts (a) and (b) are plots of the IMDS (defined by 
4>(r) = 0) of confined perpendicular Lj_ lamellae. In (a), the 
two walls favor the B monomers, a ± = 0.2hqo<j)q > 0, and the 
B domains are larger than the A domains at the walls. In (b) 
<t~ = — a + = 0.2hqo4> q , and the A domains have large size at 
y = —L, while the B domains are larger at y — L. The Flory 
parameter is N% = 11 and r s = O.lhqQ. 
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Fig. 4. Schematic drawing of the IMDS lines. The confining 
walls are at y = ±L = ±do- The dotted line is the location of 
the unperturbed IMDS. The lateral deviation from this line at 
the walls is Ax. The angle between the tangent to the IMDS 
and the a;-axis is 9. 




Fig. 5. Temperature dependence of the shape of the IMDS. In 
(a) the Flory parameter is \N = 12 (relatively strong segre- 
gation), and the IMDS are nearly flat. As the temperature is 
raised and approaches the ODT, \N = 11 in (b), the lamellae 
can easily deform in accordance with the surface fields a ± . The 
shape of decaying undulations is given by Eq. (Gq) with c = 0. 



The parameters chosen are a ± = 0.5/igo and r s = 0. 
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Fig. 6. Concentration profiles for the confined parallel Ly 

Fig. 7. Film free energy per unit area as a function of inter- 
phase. Dotted line [<f> — 4>t(y)] and solid line [4> — (j> b (y) + 5(j>(y)] 

plate separation 2L. Shown are the bulk approximation to the 

are normalized by <f> q . In (a) the two walls favor the B 

free energy of the L± phase (horizontal dashed line) , bulk ap- 

monomers (<f> < 0), a ± = 0.5hqo<f) q > and the film is sym- 

proximation of Ly (dotted line), full free energy of the L± (solid 

metric, while in (b) the film is asymmetric, a~ = — a + — 

line) and of the L\\ phase (thick solid line). The film is taken 

— 0.5hqQ(f>q, and the A monomers are adsorbed at the y — —L 

to be symmetric. In (a) a ± = QAhqf^q, while in (b) the sur- 

wall. The bulk Flory parameter is \N = 10.6 and its surface 

face interactions are smaller, a = 0.2hqo4> q - Free energies are 

modification is t s = 0.125/igi). 

measured with respect to the free energy of the bulk lamellar 
phase. Inset in (a) is an enlargement of the L± free energy by 
a factor of 10 3 , showing a deep minimum for 2L < do. The 
bulk Flory parameter is \N =11 and its surface deviation is 
t s = 0.35hq$. 
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Fig. 8. The stability of Ln (in light) vs. L± lamellae (in dark), 
as a function of wall separation 2L and interfacial strength 
a + . The free energies are taken from Eqs. (^) and ([^8|), re- 
spectively. In (a) the surface Flory parameter is the bulk one, 
t s — 0, while in (b) r s = 0.1/igo > 0. The Ly phase is pushed 
upward in the stability diagram in (b), removing the degener- 
acy between L± and L\\ that occurs for neutral walls (cr = 0) 
when t s = 0. The calculation is done by the bulk approxima- 
tion, cf>(r) — 4>b{r). The ratio a + /a~ = —2 is kept constant 
and the Flory parameter is xN = 11. 
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Fig. 9. Stability of Lm (in light) vs. L± lamellae (in dark), 
comparing in (a) the bulk approximation [free energy equa- 
tions (^) and (Bq)], with the full free energy in (b) [calculated 
from Eq. (|l|) with Eqs. @ and respectively]. Note that 
the large a + behavior lies outside the range of validity of our 
linear model. In both parts (a) and (b) t s — 0.3/w/o, the Flory 
parameter is \N = 10.8 and the film is symmetric, a~ — a + . 



